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q , We compute the triple Pomeron vertex from the Wilson line formalism, including both planar 

and non-planar contributions, and get perfect agreement with the result obtained in the Extended 
Generalized Logarithmic Approximation based on reggeon calculus. 
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I. INTRODUCTION 



It is well known that the Froissart bound [T[ for total cross-section, a tot < 7r jm\ In 2 s [2], is violated in 
perturbation theory within the Leading Logarithmic Approximation (LLA). The hard Pomeron exchange ob- 
tained in the LLA is described by the BFKL equation [3j-|6[ . This violation persists also at the Next to Leading 
Logarithmic Approximation (NLLA). Although strictly speaking valid only for hadronic observables, and not 
for external virtual states such as 7*, there is a common belief that this bound should be satisfied within any 
reasonable perturbatively resummed scheme. This was the starting point of various lines of research, which led 
to various unitarization and saturation models in which the triple Pomeron vertex is a key building block. 

The Generalized Leading Log Approximation, which takes into account any fixed number n of f-channel 
exchanged reggeons, leads to the Bartels, Jaroszewicz, Kwiecinski, Praszalowicz (BJKP) equation [7j-Hoj. The 
BJKP equation can be reformulated as a 2-dimensional quantum mechanical problem with n sites, each one 
corresponding to one of the (gluonic) t— channel reggeons (in the coordinate space), with time ~ Ins. The 
underlying hamiltonian is holomorphically separable and invariant under global conformal transformations (ill . 
H2| . In the lar ge N c limit, this 2-dimensional quantum mechanical model greatly simplifies, and turns out to be 
integrable |12h16| . The singlet color bound states of an even number n of t— channel reggeons have the quantum 
numbers of the Pomeron (with P = C =+1), while for n odd, such bound states contribute both to Pomeron and 
Odderon (with P = C =—1) exchange. For Odderon, the obtained trajectory satisfies «o < 1 When 
summing with respect to n, it is is expected that the whole series, although divergent, could have a critical 
behavior with an Odderon intercept a<Q = 1. However, these bound states decouple from Born impact factors. 
They couple to photon impact factor only through non trivial color states, at least of quadrupole type, which 
are therefore suppressed by 1/iV 2 powers. In contrast, it is possible to exhibit a critical solution (a© — 1) which 
couples to Born impact factors. These peculiar solutions can be obtained either from the perturbative Regge 
approach [23[ or from the dipole model, see Ref. [24j . 

In comparison to the previous approach, the Extended Generalized Leading Log Approximation (EGLLA) [25l — 
[29( | . in which the number of reggeon in t— channel is not conserved, satisfies full unitarity (in all sub-channel) 
and leads to an effective 2-d field theory realizing the Gribov idea of Reggeon field theory 30] in QCD (for a 
pedagogical review on this approach see Ref. [29| )■ In the framework of EGLLA, the simplest new building 
block (with singlet sub-channels) is the triple Pomeron vertex [2614281 l3l| . The conformal properties of this 
vertex allow one to relate it to the conformal blocks of an underlying (still unknown) conformal field theory, 
and using bootstrap properties, it was possible to evaluate this vertex j32|. The Pomeron vertex contains two 
contributions: a planar one, and a non-planar one, which is suppressed by a factor of 1/iV 2 with respect to the 
planar one: 

' ^planar ^non — planar ' V / 
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The dipole model [33j,|34|], [35l438j . equivalent to the BFKL equation at LLA [38, 39], is based on the description 
of the wave function of an onium state in terms of color dipole degrees of freedom in the 't Hooft limit and at 
?e s. This wave function satisfies a non linear evolution equation. In this approach the vertex of 1 2 dipoles 
-|42j is equivalent to the planar part of the triple Pomeron vertex derived through the EGLLA approach. 
The non-planar part of the triple Pomeron vertex cannot be derived from the dipole model since this approach 
relies on the large N c limit which suppresses all non-planar contributions. 

Since the time it became clear that at high-energy (Regge limit) non-linear effects dominates the dynamics 
of the scattering processes, non-linear evolution equations were derived. One of these equations is the Balitsky- 
Kovchegov (BK) equation, derived first by Balitsky (43l - l46| in the Wilson line formalism, and then by Kovchegov 
[H SI] in the dipole model. 

The Wilson line formalism is an operator language. It is based on the concept of factorization of the scat- 
tering amplitude in rapidity space and on the extension to high-energy (Regge limit) of the Operator Product 
Expansion (OPE) technique, which was known before only at moderate energy (Bjorken limit) as an expansion 
in terms of local operators or in terms of light-ray operators. In Deep Inelastic Scattering (DIS) off a hadron 
at high-energy, the matrix elements made of Wilson line operators appearing in the OPE, describe the non 
perturbative part of the process, and their evolution in rapidity is related to the evolution of the structure 
function of the target. In order to find the evolution equation one may use the background field technique. 
The Wilson-line operators evolve with respect to rapidity according to the Balitsky equation, which reduces 
to the BK equation in the large N c limit. The BK equation describes the so called fan diagrams neglecting 
all non-planar contributions, while the Balitsky equation gives one the possibility to describe the fan diagrams 
including also the non-planar contributions. 

A similar method to the background field technique is given in Ref. [49], where the authors studied the 
propagation of a fast moving particle, and showed that resumming the emission of soft gluons from this source 
and neglecting its recoil (eikonal approximation), one may obtain the BFKL equation. 

The Color Glass Condensate (CGC) [50T458T ] (for a review, see Ref. [59]) is another available way to describe 
high-energy scattering processes. This approach is similar to the Wilson line formalism described above, and 
the corresponding evolution equation is the Jalilian-Marian, Iancu, MacLerran, Weigert, Leonidov, and Kovner 
(JIMWLK) equation. Indeed, the Balitsky equation, as already mentioned above, is not a closed evolution 
equation, but it contains an infinite set of evolution equations which goes by the name of Balitsky-hierarchy. 
The JIMWLK equation is equivalent to this hierarchy of evolution equations thus leading to the acronym 
B- JIMWLK equation. 

In all of these approaches, the description of the scattering of two probes involves the computation of the 
interaction of one Wilson loop describing one probe with the field of the other. Note that within the dipole 
model, a similar approach, based on the computation of the scattering phase of a dipole in the field emmited 
by a fast moving object (involving color structures as well as multicolor states), was suggested in Ref. [60]. 

A further a ppr oach to the description of high-energy processes is based on an effective field theory. This 
formalism [6l]464j provides the building blocks (reggeon-reggeon-gluon) necessary to the explicit computation 
of any type of diagram at this regime. However, the precise relationship between this effective theory and the 
EGLLA approach has not been clarified yet, and explicit applications of this effective field theory is highly 
desirable [65146 71 ] . 

The triple Pomeron vertex is the first non-trivial building block common to all the above approaches. It turns 
out that up to now, its exact expression, including planar and non-planar contribution, was derived only in the 
EGGLA approach 28]. Using its conformal invariance |3l|. both the planar and the non- planar contributions 
of the vertex were computed in the coordinate space [68| (see Ref. [32j for explicit expressions in the SU(N C ) 
case). 

So far, only the planar contribution was derived with formalisms different from the EGLLA one. These are 
the dipole model [4014421 ] and the Wilson line formalism applied to diffractive processes [o9j . 

The purpose of the present paper is to show that the exact expression of the 3 Pomeron vertex (planar and 
non-planar contributions) can be easily derived through the Wilson line formalism not only for diffractive case 
but also for fan diagrams. 

Indeed, this formalism was already used to derive several new and desirable results, and to confirm others 
which have been derived after many years of calculations. This include for example the NLO BK kernel in 
QCD [7(3 and in Af=4 SYM [jlj whose linearized version confirmed the NLO BFKL kernel in QCD [72411 



and in A/"=4 SYM [76J. Moreover, through Wilson line formalism it was possible to derive for the first time the 
conformal expression for the NLO BFKL kernel [77j, the full amplitude at NLO in A/"=4 SYM for four scalar 
currents made of chiral-primary operator [78|, and the analytic result for the NLO photon impact factor [79j . 
relevant for phenomenology at high-energy, was recently derived. 
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It is clear that the Wilson line approach to the study of high-energy scattering processes opens the way 
to attack more difficult problems which have not been solved before despite the many efforts that have been 
devoted to them using different techniques. This includes for example multiple Pomeron or Odderon vertices, 
as well as sub-leading contributions. 



II. INTRODUCTION TO THE WILSON LINE FORMALISM 



In this section we will give a brief introduction to the Wilson line formalism. This will set the notations 
which will be used for the derivation of the triple Pomeron vertex first in the diffraction case and then for fan 
diagrams which are important for unitarization. 

The main tool which we will use is the OPE for high energies [43[ of the T— product of two electromagnetic 
currents in terms of Wilson lines: 



Tj fl {x)j u (y)= / d 2 z 1 d 2 z 2 I^(x,y;z 1± ,z 21 _)Tr{U{z 1± )U^z 21 _)} + --- (2) 



with the operator j M (a;) = ^(x^^lx) . This expansion is in terms of a coefficient function to be identified with 
the photon impact factor and a matrix element of two Wilson line operators, with 

U(x x ,y x ) = l-±Tr{U(x ± )U\y ± )}, (3) 

where the Wilson line is defined as usual by the operator 

U(x±) = PexpUg [ dup^A^pxu + xM. (4) 

J — oo 

We use here the standard Sudakov decomposition k — ap% + (3 p 2 + k± , where pi and p 2 are light-like vectors 
defined in such a way that in a typical high-energy scattering process, the first projectile (refered as the "above" 
one) flies almost along p\ while the second one (resp. "below") flies almost along p 2 . 
The operator IA evolves according to the Balitsky equation 



d -ft, \ a sN c f, 2 f> - y)\ 



— U(x±,y±) = 3 _ 2 C / d 2 z — — 2 2 [U(x±,z±) + U(y±,z±) -U(x±,y±) - U(x±, z±)U(z±,y±)] ,(5) 



where x± y±, z± are two-dimensional vector with Euclidean metric 1 . The r\ dependence of the operator IA 
enters as a regulator of the divergence by changing the slope of the Wilson line according to p\ — )• p\ + e~ 2r, p 2 
in Eq. flU). 

The BK equation [43l . |47| is obtained from Eq. ([5]) at large N c when the correlation function of the non-linear 
term IA (x, z) Li(z, y) decouple in a product of two correlation functions (U(x, z)) (U(z, y)) . From now on we will 
use the short-hand notation U{x, y) = lA xy ■ 



III. TRIPLE POMERON VERTEX FROM DIFFRACTION 



A. Diffraction within Keldysh formalism 



We will derive in this section the triple Pomeron vertex for diffractivc processes, using Keldysh formalism 
adapted to describe diffractive processes through functional integration [80] - The idea is to use the OPE for 
diffractive high-energy processes [69| in order to reproduce automatically the Cutkosky rules for the calculation 
of total cross-sections. One introduces two different fields, each of them living on one side of the cut, which 
results in three different propagators: (A + A + ) , (A + A~) , (A~A~). 



Hereafter, these variables will be denoted simply like x, y, z. 



Following [69|, (8l|, we write the diffractive amplitude of the 7*p — > p r + X process as 

W diS = J2 e ^J ^^I A (kuO)(N\Tr{W i=m2/s {ki_)W u=m ^ s (-kA_)}\N}, 



flavors 



where W(k±) is the Fourier transform of 

W(x ± ) = V^(x ± )U(x ± ), W\x ± ) = U\x x )V{x x ) , 



(0) 



(7) 



and U(x±) now denotes the Wilson-line operator constructed from A + fields while V(x±) denotes the same 
operator constructed from A~ fields: 



U(x±) = Pexplig du p^ A+(p 1 u + x±)\ 

^ J — oo ^ 



V{x±) = Pexp 



du A] 



(piu + xx)\ . 



(8) 



In Eq. (J6j> m is the mass of the proton, and the notation £ = e v is used. In the case of diffractive processes, 
the operator 



J * r 



(9) 



evolves according to the same Balitsky equation J5]), as [6 



d tAii \ asNc I J* 



[x - z)\{z - yY 



{ W(x, z) + W(z, y) - W(x, y) - W(x, z)W{z, y)\ . (10) 



In the RHS of this equation, the non-linear term should be interpreted as the splitting of a diffractive Pomeron 
defined by (VV(x±, y±)) . We now define V in a similar way as U in Eq.Q 



V(x u y x ) = I - l-Tr{V(xx)VKy±)} ■ 



(11) 



Our goal is now to extract from the non-linear part in Eq. (fTUj) terms of the type (U) (V) which will be 
interpreted respectively as the Pomeron on the left (resp. right) of the cut. In order to do this, we should 
linearize this non-linear term up to two gluons accuracy 2 . 



B. Linearization of non-linear term 



The idea of linearization consists in expanding the non linear term of Eq. (1101) up to g 4 and rewriting this 
result in terms of products of the type 1A V . The two gluon approximation means that each IA and V should 
be approximated up to g 2 , since each Pomeron is a singlet color object. The obtained result, whose detailled 
derivation is given in the Appendix, is: 



Tv{W x W^ - l}Tr{W z W$ - 1} ^ U xz V zy +U yz V zx 



1 



N? 



■7 ^zx Vzy\ 



I r / 

UxzV zy ~f~ MyzVzx H~ ~j^2 Uxy^xy ~f" ^xz^xz ~\~ ^zyVzy ~r~ t4xy\^xy vxz ~ ^zy 



N*-l 

~^xy (^xy t^xz ^zy 



(12) 



2 In the following, the notation * on operators will be removed for simplification. 
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C. Projection on BFKL Green functions 

In order to extract the triple Pomeron vertex, one has to define precisely how to factorize out each of the 
three Pomeron Green function from the three- Pomeron correlator. This is achieved according to Fig. [TJ We 
denote the above (below) Pomeron Green functions by (resp. VP). The vertex y lp ^ 2P i s defined symbolically 
in the following way: 

(PPP) = (AA^')^ 1P " >2P *^ (13) 

where AA^P' is the amputated Pomeron Green function which is denoted as y p in Fig. [TJ This prescription 
is in accordance with standard definitions, as exhibited for example in Ref. [82j |. from which notations of Fig. 
[JJare inspired. We will now translate this definition of the triple Pomeron vertex in the shock- wave approach. 

In the present treatment, we deal with colorless probes. These probes are dipoles, which respect the global 
conformal invariance of the BFKL equation. The dipole-dipole scattering, in the BFKL approximation, can 
then be presented as an elementary function of a conformal anharmonic ratio. This is the basis of the so-called 
Moebius representation of BFKL. 

In the shock-wave analysis, (I4 xy ) describes the two-gluon non-amputated amplitude in the Moebius represen- 
tation (the upper probe is a dipole with coordinate x and y referring respectively to the position of the quark 
and anti-quark pair), where the average is on the external field of a lower probe. The contact with the usual 
4-gluon BFKL Green function in the Moebius representation can be made if one considers the specific case of a 
lower probe made of two Wilson lines in a color singlet state (a dipole), each of them having definite transverse 
coordinates. 

The dynamics of the process is encoded in the dipole kernel (identical to the BK kernel) which acts on the 
coordinates x and y (the same situation appears for the evolution of multiple dipole densities in the dipole 
model, which is the starting point of the Kovchegov approach) , while the averaging from below on a given probe 
does not affects this dynamics. In both the BFKL and the shock-wave pictures, this dynamics is encoded in a 
kernel which acts on non-amputated functions (this remains true also in the dipole picture). 

Now, on one side, the BFKL Green function is non-amputated, both from below and from above. Indeed, 
at Born order it simply reduces to the product of two propagators. Its Moebius representation is obtained by 
Fourier transform and then by a "substraction", which is needed to enforce its vanishing for equal upper or 
lower coordinates 82]. This can be obtained directl y w hen computing the elementary dipole-dipole scattering 
amplitude in the two gluon exchange approximation (39(). On the other side, in both the shock- wave and dipole 
formalisms, the operator U or the dipole densities correspond to amputated quantities from the point of view 
of the below probe. 

Indeed, to get a scattering amplitude one should convolute these amputated Green function to the below 
probe, thus restoring the gluonic propagators, as it is indeed done in these two formalisms. In the Wilson 
line formalism the fields are not contracted with the probe from below (no propagators). In the dipole model, 
when computing the scattering of two onia, in the frame where the below onium is almost at rest, one first 
evaluates the dipole content of the above probe. Then, the restoration of these gluonic propagators appears 
when convoluting this dipole density with the elementary dipole-dipole scattering amplitude between an internal 
dipole constituent of the above probe and the below dipole. 

D. Projection on conformal three-point functions 

In the Moebius representation the Hamiltonian of the BFKL equation acts as an integral operator, according 

to 

tt H BFKL f{x, y)= J d 2 z [x } X z ~^ 2 _ z)2 [f{x, y) f{x, z) fiz, y)] . (14) 

Based on the conformal invariance of this BFKL hamiltonian, one can diagonalize it using the following set of 
eigenfunctions 

^(Pio,P 20 )^(-^) h (^f) h , (15) 
VP10P20/ VP10P20/ 



Figure 1: A typical diagram contributing to the triple Pomeron vertex V, with three 4-gluon Green functions (denoted 
as ^> and W in the Moebius representation) at g 2 order. This exhibits explicitly the amputation of the above Green 
function from below. 





b) 



Figure 2: a) Configuration of planar diagrams, b) Configuration of non-planar diagrams. 



where the conformal weights h and h are given by 

, 1 + n . 
h = — h i v , 



— 1 — TL 

h = 1 - h* = — Yiv. 



(16) 



The eigenvalue equation then reads 



with a 8 = as ^ c , and where u(h) is given by 

u(h) =2a s Re[ip(l) -ip(h)} . 
The corresponding BFKL Green function then reads, in the Mellin space w conjugated to s, 



G%™*{x, y, x', y') = E 



h{h-l) 



(17) 



(18) 



d 2 Po E h - h {x -po,y-po) KAx' - Po ,y' - p ) ■ (19) 
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Let us now turn to the shock-wave formalism. In the BFKL approximation, the average of the 11 operator on 
a given probe (denoted as B) can be written as 



; B| U%> |B) = 5>. y | a) (a | B) L— 



h a (h a - 1) 



(20) 



where the notation \a) is used in order to label both the center of mass p a and the conformal weights (ft a , h a ). 
Note that in Eq. (f2"0"f we make an identification of the evaluation of the operator U on the state B through 
functional integration on the LHS, with quantum mecchanical notations appearing on the RHS. Here \x,y) 
denotes an upper dipole, x and y being the coordinates of the two corresponding Wilson line. The relationship 
with the 4-reggeon BFKL Green function G BFKL (x, y, x', y') in the Moebius representation is obtained when 
choosing the lower state B to be a dipole of coordinates (x', y') , for which (TK)1) reads, in the above notations, 



( dipole (x', y') | U<$ | dipole (x', y') ) = E ( x > y\^){a\x',y') 



h - i 



2 



h a (h a - 1) 



(21) 



w — uj{h a ) 

leading to the following natural identification: 

E h a ,h a ( x ~ P<*> y-Pa) = (x,y\a), E* a fi a (x'-p a ,y'-p a ) = {a\x',y'). (22) 

We denote the quantum numbers of the upper Pomeron by a , while the lower left Pomeron (resp. lower 
right) is labelled by j3 (resp. 7) (see Fig. We are interested here in the extraction of the triple Pomeron 
vertex. For that purpose, based on the factorised form ([21]) of the BFKL Green function as a series of products 
of conformal blocks, the two lower Pomerons are described through (x,z\fi) and (2, y | 7 ) while the above one 
is described through (a | x, y) . Now, using the property that 

{E h ,h{x ~ P,y ~ p))* =E 1 _ hA _ Jl (x-p,y-p), (23) 

we will not make any distinction between the above and the two below Pomerons, which will all be described by 
a set of E h ^ functions, as it is done for example in Ref.[32l| to which we would like to compare our final results. 

We now fix our notation for the internal degrees of freedom of the vertex, based on the non-linear term of 
Eq. (I10|) . Let us equate the coordinates x, y and z with p a , pb and p c respectively. It is convenient to identify 
the component of the dipole with definite conformal weight with U , according to 

U xy -> E hft f te (p a p,p b0 ), U xz -> E hp j lfj (p a0 ,p c p) 

V xy -> E h ^ (p ai , p^), V yz -> E hi - hi (p b7 ,Pcj) ■ (24) 

Now, since we are in the Moebius representation, for each of the three Pomeron which are here involved, the 
Green function is the sum of products of conformal block with even conformal spin n. This is due to the fact 
that 

G^ b L ius (x, y, x', y') oc G BFKL (x, y, x> , y>) + G BFKL (y, x, x' , y') , (25) 

as discussed in Ref. [39| (see Eqs. (2.50) and (2.51), as well as Eq. (B.12) in the Born approximation). Therefore, 
since n = h + h is even, and using the property 

E h - h (x -p,y-p) = (-l) h + h E hrh (y -p,x-p), (26) 
one can make the following identification 

U a b = Uba = E h h (Paa, Pba) = E h„h n (pben Paa) (27) 



where ti a b = Uba in the two-gluon approximation. Identifications similar to Eq. (|27l) are done also for £4c> Vab 
etc... 

Since the upper Pomeron is amputated, one needs to evaluate the effect of the amputation on a conformal 
block. This is obtained using the following identity: 

/^ a /^ b E ha t la {p aa , p ba ) = 16h a (h a - l)h a (h a - l)E ha i a (p aa , Pb a )~2~ =2~ • ( 28 ) 

Pab Pab 
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The three-pomeron vertex is now denned, according to Eq. (|12p , as 
1 



>2P 



d pad p b d p c 16 h a (h a - l)h a (h a - l)E h % (p aa , p ba ) 



\Pab\ 



PlbPlb \Pac\ 2 \pbc\ 2 



UacVcb + UbcVc 



1 



-u. 



'abVab + MacVac + UcbVcb + Uab (V ab - V ac - V cb ) +Vab{U ab 



~ M ac ~ U c b 



Using property (j2"T|) , the planar part in Eq. (f2"!))) can be rewritten as 



V, 



1P^2P 



planar 

16h a (h a - l)h a (h a - 1) / d 2 p a d 2 pbd 2 p c E h z (pba, Paa) 



(29) 



(30) 



\Pab\ 2 \pac\ 2 \pbc\- 



-, E hf,hf,{PaP: Pcp)E h ^(p cl , Pb-y) ■ 



Let us now consider the non-planar part of the vertex. Using the hermiticity property of H BFKL (TT4"]) when 
acting on the product of the two below Green functions U xy V X y and acting on the above amputated Green 
function we have, 

7T J d 2 Pa d 2 pbd 2 p c E ha - h J Paa ,p ba ) [ff™ U a bV a b] 

= tt / d 2 p a d 2 p b d 2 p c [H^ b FKL E h j la {p aa , p ba )]E h - h (p a0 , p b p)E h z (p ai , p bl ) 



= 7T / d Pad p b E haha (p aa , Pba) Eh B h 8 (pa/9, Pbp)E h ^ (/J a7 , Pfe 7 ) 

Similarly, for the other two terms of the non-planar part we have 



-u)(h a ) 



(31) 



d p c - 



I Pah | 



Pac\ 2 \Pbc\ 2 



Uab^ab - Vac ~ V cb ) + V ab (U ab ~ U a c ~ U c b)] = 7T U a b [H^ FKL V a b] + ^ V ah [fff b FA ' L W ab ] 



E hghg{pa0,Pbp)H a b E h^h-,(Paj,Pb-y) + E h % (fiaj, Pb~/)H ab E h h {p a p , p bj j) 



= -TT— [w(hp) +uj(h^)]E hff f lfj (p a0 ,p b& )E h ^(p a7 ,p bj ) 
Thus, using once more property f[27|) , the non-planar part reads 



(32) 



V, 



1P-S.2P 



2tt 



non-planar ~ ~ Jf2 J d Pa d Pb lQh a (h a - l)h a (h a - 1 )E ha ha (pba , Paa) 

x r— jjRe^Cl) + i>{h a ) - t/j{h/3) - ip(h^)}E h?Jl0 (pbfs, p a p)E h % (p^, p ai ) . 

\Pab\ 



(33) 



Finally, putting together the non-planar part Eq. ([3"3"|) and the planar one Eq. (j3"0"l) . the triple Pomeron vertex is 

rlP^2P 

I 



V 

d 2 p a d 2 p b 16 h a (h a - l)h a (h a - l)E h t (p ba , p aa ) 



d p c 



2tt 1 



AT2 I 14 M^ 1 ) + V'('la) ~ 1p(hp) ~ 1p(h 7 )}E h (p b p,p af} )E h ^ (fty, p^) 

^c \Pab\ 



(34) 



We now compare our result (|3"4")l with the one obtained through the reggeon approach of Ref. [28[ , and which 
is written explicitly for SU(N C ) and for arbitrary conformal weights in Ref. [32j. Up to a global normalization 
factor, related to the convention used to define the triple Pomeron vertex, Eq. (f3"4")l agrees with Eq. (2.1) of 
Ref. dl. 
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IV. FAN DIAGRAM APPROACH 



The diffractive case, discussed in the previous section, implied that the upper Pomeron was at t — 0. In this 
section instead, we will show how to obtain the triple pomeron vertex without using the Keldysh formalism. 
This will allow us to obtain the triple Pomeron vertex for fan diagrams which means that every Pomeron is now 
at arbitrary t. Therefore, we need to work with the Balitsky equation ([5]) which we rewrite here for convenience 

-j- U{x,y) = / d 2 z ^ [U(x,z] + U(y,z) -U(x,y) -ti(x,z)U(z,y)] . 

dn 2tH J [x- z)j_{z - y) z ± 

It is easy to see that at large N c limit the correlation function ( U(x±, z±) U(z±, y±)) decouples to the product 
of two correlation functions (hi(xx, z±)) (Vi(zx,y±)) and this non-linear term is interpreted as the splitting of 
one Pomeron into two Pomerons. In this way the Balitsky equation with the truncation of the hierarchy reduces 
to the BK equation; its non linear term coincides exactly with the planar part of the triple Pomeron vertex [82j . 
The triple Pomeron vertex takes the following form 

V™ cx / d*z { J X z -J^_ y)l (U xz )(U zy ) . (35) 

Our aim is now to extract from the non-linear term U xz lA zy , not only the planar contribution to the triple 
Pomeron vertex but also the non-planar one. We then show that the result so obtained coincides with the 
one we obtained for diffractive processes. We will adopt the following procedure. First, we consider the 
correlation function of four Wilson lines i.e. (L( xz U zy ), we then apply the two gluon approximation to them 
and finally we rewrite the contributions thus obtained in terms of decoupled correlation function of the type 
{U xz )(U zy ) . This method is technically very similar to the one used in Sec. IIII1 the details of which are given 
in the Appendix. We first expand each Wilson line operator. In what follows, we use the shorthand notation 

U x k, U x ^ + Ux + U x 2 ^ + . . . , with being the zeroth order term of the expansion, Ux the first order term, 
and so on. Thus, we have 

N 2 c (U xz U zy ) (36) 
« (Tr{l - (t/J°) + E/« + Uj»)(UW + f/«t + U^)} Tr{l - (E/J°) + £/« + U^)(U^ + + U™)}) . 

At this point we want to rewrite eq. (|36| as a sum of decoupled correlation functions like {U xz )(U zy ) such 
that when we apply the 2 — gluon approximation to them we get back the Eq. (|36[) . Since E/(°) = 1, the only 
possible contraction in order to produce terms of the type {U xz ){U zy ) are between terms or between 
terms. Contraction between terms of order higher than 2 would result in remaining multiplicative terms which 
are not color singlet (these terms are of the type TrU^ which vanish). Other terms involving contraction of 
gluon fields at the same coordinate will clearly not contribute: terms like U zz vanish. Finally, the expansion of 
each Wilson line is needed only up to second order. So, it is then easy to see that 

(U xz U zy ) (37) 

^ 2 (J-l) { 2{Uxz){Uzy) + l^[ 2{Uxy) ( {Uxy) ~ {Uxz) ~ {Uyz) ^ + { U "»K U *») + (^*>< W ~) - (U*v)(U*y) 

We can immediately recognize in Eq. (|38[) the planar contribution {Li xz ){U zy ) which coincides with the non- 
linear term in the BK equation and with the planar part of the diffractive triple Pomeron vertex obtained in 
the previous section (cf. Eq. (|47p). The terms proportional to N~ 2 are instead the non-planar contributions 
which are suppressed in the limit of iVc — >• oo and which are, therefore, obtained only from Balitsky equation 
and not from BK equation. 

We now want to compare the two above approaches: the one based on diffractive processes, Eq. (|29|) . and 
the other one based on the fan diagram approach, Eq. (|38p . The first obtained result is a particular case of the 
second one, since it was derived for the splitting Pomeron at t = . Indeed, the second one can be obtained when 
identifying hi with V since in the fan diagram case one cannot distinguish between the two produced Pomerons. 
Let us first consider the planar contribution. In the diffractive approach, the obtained structure (see Eq. (fT2"jl ) 

N 2 

j^j[UacV cb + V ac U cb ] (38) 
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is an operator which should be contracted with an external set of Pomeron states of quantum numbers denoted 
by \(3) and I7) (which as already stated above describes both the conformal weight and the center of mass 
coordinate of the Pomeron state). At this stage, in the case of the diffractive amplitude, these two Pomeron 
states are distinguishable (one is at the left of the cut while the other one is at the right), and after using the 
symmetry of the integrand under the replacement a <-» b, the net result reads symbolically 

' X " " l/?)V cb | 7 ). (39) 



N*-l 

Now, one can make the identification of the £4/|/3) and Vij\/3) states, leading to the final result 

N 2 

~ ' Z4 c |/3)Z4 fc | 7 }. (40) 



m - 1 



On the other hand, from the fan diagram approach, one has (see first term of Eq. (j38|) ). using the same overall 
normalization as in Eq. (|38[) . 

iV 2 N 2 

■2W ac W cb |/3)| 7 ) = — ^2(U ac \l3))(U cb \j)) (41) 



2(N 2 - 1) "^"w "" JVf -1' 



where the factor of 2 on the RHS of Eq. (|41|) is due to the two possible contractions. This shows explicitly that 
both, the planar result obtained from the general fan diagram case and the one obtained from a continuation 
of the t — diffractive case in Keldysh formalism, are in agreement. The proof for the non-planar case follows 
the same line of thinking. 



V. CONCLUSION 



In this paper we have shown that the triple Pomeron vertex, including the planar and the non-planar con- 
tribution, can be obtained very easily within Wilson line formalism. In Ref. [691 ] this was already done for the 
case of diffractive processes using Keldysh formalism, but the result obtained there was only for the planar part 
of the vertex. In section III, we have shown how to compute also the non-planar contribution of the diffractive 
triple Pomeron vertex from Wilson line formalism. To this end we considered the generalization of the Balitsky 
equation for diffractive process, and using the linearization procedure and the 2-gluon approximation, we have 
extracted the desired subleading term in N c of the vertex. 

In section IV we have extended the result of section III to the more generic case of fan diagrams, where the 
Pomeron, which split to two other ones, does not need to be at t = 0. We then showed that the triple pomeron 
vertex for fan diagrams is the same as the one obtained in the diffractive case. 

Since, as we have shown in the present paper, the Wilson line formalism allows one to re-derive very easily 
results that have been obtained after non trivial and lengthy calculation, we plan to use it to study other, still 
unknown and highly desirable results. For example, it will be interesting to compute the vertex for P — > 3P [83|, 
P->00 [84], O — s- PO (so far unknown), or more generally nP — > mP (inaccessible at the moment through 
reggeon calculus techniques). These non trivial building blocks will be relevant in order to identify the unknown 
underlying effective theory for high-energy scattering processes. This study is in progress. 
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Appendix 

Below we present some details of the linearization procedure which we used to extract the planar and non- 
planar part of the triple Pomeron vertex. 

Let us consider the non-linear term in Eq. (fTOj) 

N 2 Tr{W(x, z)W(z, y)} = Tr{W x W f z - l}Tr{W z W} - 1} = Tr{V z V^U x U f z - l}Tr{V y V}U z U y ~ 1} ■ (42) 

Our aim is to extract from this expression the contribution of two non-interacting Pomerons: one Pomeron 
built by v4 + -fields and the other one built by ^"-fields. We now approximate Wilson line operators U and V 
up to linear terms in A + and A~ fields (here we use the short-hand notation iA+ = ig J Q dup^A^ (piu + x±) 
and similarly for iA~) 

N 2 Tr {W(x,z)W(z,y)} 

« Tr{(l + iA~)(l ~ iA-)(l + iA+)(l - iAt) - l}Tr {(1 + iA y )(l - iA~)(l + - iA+) - 1} . (43) 

From each trace of Eq. f|43[) we keep only terms to second order since we are working in the 2-gluon approx- 
imation. Thus, we have 

N 2 Tr{W(x,z)W(z,y)} « Tv{A^ A~ + A+A+ - A~A+ + A~A+ + A~ A+ - A~A+} 

x Tr{A-A- - A- At + A~A+ + A z At - A Z A+ - At A+} . (44) 

The next step of the linearisation procedure up to 2 gluon accuracy consists in keeping in the product of the 
two traces in (l4"4l only terms which involve two A + fields with different coordinates and two A~ fields with 
different coordinates. Thus, we obtain that 

N 2 T r {W(x,z)W(z,y)} % Tr{A-A;}Tr{A+A+} + Tr{A+A+}tr{^;Aj} + Tr{A;A+}Tr{A;A+} 
-Tr{A-A+}Tr{A~A+} + Tr{A; A+}ti{A; A+] Tr{A^ A+}Tr{A; A+} + Tr{A~ A+}Tr{A- A+} 
+Tr{A-A+}tx{A-A+} - Tr{A~ A+}Tr{A~ A+} - Tr{A^A+}Tr{A^A+} + Tr {A~ At }tr{A~ A+} . (45) 

We now have to rewrite Eq. (|45[) in terms of product of traces involving only A + and A~ fields, so we obtain 
N 2 Tr{W(x, z)W(z, y)} % Tv{A-A:}Tr{A+A+} + Ty{ At At }Tr{ Ay A z } 

+ [Tr{ ^ A y } ~ Tl ^ A - A ^ Tr i A t A t}] [T*i A *^} ^{ A * A z} Tr{A-A;}] . (46) 

The final step is to write Eq. (|46| in terms of the original operators U and V such that when we apply to them 
the 2-gluon approximation we get back Eq. (j46]). So, we have 

N 2 C Tr{W(x, z)W(z, y)} « N 2 (u xz V zy + U yz V zx + j^—j [U xy - U xz - U yz ] [V xy - V zx - V zy ]j (47) 

which is the linearization we used in Eq. (|12[) . 
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